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I. Introduction et Justification du choix : 


Lors de notre cours de béton armé on a pu savoir la nécessité d’utiliser des 
armatures : 


Le béton possède une grande résistance à la compression et une résistance 
moindre à la traction. Dans les structures en béton se développe un ensemble de 
contraintes générées par les diverses actions auxquelles elles sont soumises. La 
résistance à la compression du béton lui permet d’équilibrer correctement les 
contraintes de compression. En revanche , du fait de la relative faiblesse de sa 
résistance à la traction, il n’en est pas de même pour les contraintes de traction. 
C’est pourquoi l’on dispose dans les parties tendues d’une pièce en béton, des 
armatures (barres ou treillis soudés) en acier (matériau qui présente une bonne 
résistance à la traction). Chaque constituant Joue ainsi son rôle au mieux de ses 
performances : le béton travaille en compression et l’acier en traction. 


La quantité d’armatures et leur disposition, dictées par la répartition des 
contraintes, résultent de calculs qui font appel aux lois de comportement des 
matériaux. Les bétons sont en majorité employés en association avec des 
armatures en acier. La quantité d’armatures et leur disposition, dictées par la 
répartition des contraintes, résultent de calculs qui font appel aux lois de 
comportement des matériaux. Les bétons sont en majorité employés en 
association avec des armatures en acier. 


Ainsi que les Principes du calcul du béton armé : 


Les règles de calcul sont conçues de façon à garantir la sécurité et la pérennité 
des structures. Ils précisent le niveau maximal des actions pouvant s’exercer sur 
un ouvrage pendant sa durée d’utilisation. 


Pourtant, Un bâtiment en béton armé est dimensionné à l’explosion accidentelle 
du matériau pyrotechnique qu’il abrite. Cette explosion génère une onde de 
souffle qui se propage à l’intérieur du bâtiment. Les parois d’about qui sont 
uniquement constituées de bardages sont soufflées. La structure en béton armé 
peut être endommagée mais ne doit pas s’effondrer suite à cette explosion. 


Le but c’est de suggérer un dimensionnement de ferraillage avec un cout minime 
ce qui est l’objectif de toute optimisation. 






































Il. Eclaircissement de la problématique : 


Lors de la conception d’un bâtiment pour résister au souffle résultant d’une 
explosion, la littérature normative permet l’utilisation de méthodes de calcul 
avancées. L'utilisation de méthodes normatives basées sur une approche globale 
et statique équivalente conduit généralement à un dimensionnement très 
conservateur. Une méthode de calcul avancée discrétise le phénomène 
dynamique en pas de temps (au sens numérique et physique) et la géométrie de 
la structure en plusieurs zones. Cette approche permet de restreindre les zones 
de ferraillage très dense, alors qu'une approche statique équivalente nécessiterait 
le même renforcement dans des zones plus étendues. Parallèlement, le calcul 
dynamique transitoire permet de représenter finement les phénomènes locaux de 
déformation des panneaux des éléments structuraux et de torsion des raidisseurs 
grâce à l’application d’un champ de pression variant dans le temps et dans 
l’espace. 

L’onde de surpression se propage dans le local depuis sa source, et est réfléchie 
par les parois et le plafond. Les ondes réfléchies se recombinent entre elles et 
génèrent finalement sur les parois et le plafond une succession de phases de 
chargements positives et négatives. De plus, comme la source a une géométrie 
rectangulaire et est de faibles dimensions par rapport au bâtiment, le front 
d’onde initial est ellipsoïdal. Par conséquent, même si le profil temporel de 
pression est uniforme sur l’ensemble du front d’onde, les chargements des parois 
et de la toiture sont donc des charges roulantes. La combinaison de l’ensemble 
de ces facteurs implique un chargement temporel variant selon les zones d’une 
même paroi. Nous verrons plus loin dans cet article que cet aspect est important 
pour analyser les réponses globale et locales du bâtiment. 

L’Eurocode 2 dicte les règles pour l’analyse structurale et pour la vérification 
des états limites ultimes de résistance [1]. Les phénomènes dynamiques ne sont 
pas abordés spécifiquement. La mise en pratique de l’Eurocode 2 par une 
approche conventionnelle nécessite toutefois de traduire le chargement 
dynamique en charge statique équivalente [2]. Compte tenu de la multiplicité 
des zones de chargement, des différences probables de fréquences propres selon 
les zones, et de la complexité des profils temporels de pression, cette démarche 
peut d’une part s’avérer très conservatrice, et d’autre part manquer de déceler 
des phénomènes locaux nécessitant des dispositions constructives particulières 
comme la torsion locale par exemple. 


IT. La solution proposée : 


Nous proposons une approche de dimensionnement qui permet d’appliquer les 
règles de l’Eurocode 2, tout en discrétisant l’espace et le temps. Il s’agit de 
coupler un calcul dynamique aux éléments finis linéaire avec un calcul de 
ferraillage Eurocode 2 sur le torseur enveloppe pour chacune des mailles du 
modèle. Cette méthode est applicable à l’ensemble des zones du bâtiment, y 






































compris celles qui présentent des singularités géométriques. Sur les zones 

courantes du bâtiment, le ferraillage peut encore être optimisé en appliquant 

une méthode de calcul non linéaire, dans laquelle le béton et le ferraillage sont 

explicitement modélisés avec les non-linéarités de leurs lois de comportement. 
L’aspect novateur de la méthodologie proposée est l’enchaînement du calcul 
dynamique de dimensionnement, qui permet d’obtenir une première définition 
des sections d’armatures à mettre en place, avec un deuxième calcul plus fin 
pour optimiser ces sections dans les zones où cela présente un intérêt. On 
optimise ainsi le dimensionnement tout en évitant de réaliser plusieurs fois un 
calcul fin sur l’ensemble de l’ouvrage, ce qui serait plus lourd à mettre en 
œuvre. 


IV. Explication de la méthode : 


Le bâtiment qui servira d'illustration pour l’application de la méthode proposée 
est un hall en béton armé constitué cinq portiques. Il a une longueur d’environ 
40 m, une largeur d’environ 20 m et une hauteur d’environ 15 m. La trame des 
portiques est d’environ 6 m. Il est composé de voiles et d’une dalle de toiture 
d’épaisseur 40 cm, ainsi que de poteaux de renforts de section 50 cm x 90 cm 
et de poutres de section 50 cm x 160 cm. Les poutres ont une retombée sous 
plafond de 40 cm de hauteur . 


Les voiles présentent des ouvertures correspondant à des portes. Un petit 
bâtiment que nous nommerons « chicane » est mitoyen au voile Est. La face 
sud est constituée d’un bardage considéré comme étant soufflé lors de 
l'explosion. 
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e Méthode linéaire : 


Les états de contraintes sont extraits sur chaque maille du modèle, à chaque pas 
de temps défini au préalable. Ils sont convertis par le logiciel STFER© en 
tenseurs (Nxx, Nyy, Nxy, Mxx, Myy, Mxy, Qx, Qy) [5]. 

Les armatures longitudinales sont calculées à l’aide de la méthode de Capra et 
Maury [6] dont le principe est de projeter les efforts des éléments de coques sur 
des facettes perpendiculaires à l’élément et d’effectuer pour chaque facette un 
calcul de béton armé en flexion composée. Chaque facette est définie par sa 
normale qui fait un angle Q avec l’axe local x de l’élément, et est soumise à un 
effort normal NQ et un moment de flexion MQ. Un calcul en flexion composée 
est effectué pour obtenir les forces Fb(Q) et Ft(Q) qui doivent être repris 
respectivement par les armatures inférieures et supérieures : 

FQ,b < Rb(Q) et FQ,t < Rt(Q) (1) 

L’équilibre de toutes les facettes implique l’équilibre de la coque. L’effort 
normal et le moment de flexion sont obtenus par la rotation du tenseur des 
efforts et moments : 

NQ = Nxxcos2Q + Nyysin?2Q + 2NxycosQsinQ (2) 

MQ = Mxxcos?2Q + Myysin2Q + 2MxycosQsinQ (3) 

De la même manière, l’effort résistant des lits d’armatures placés respectivement 
dans les directions x et y est obtenu par : 

AXsxcos?Q et AYsysin2Q (4) 
où AX et AY sont les sections de ferraillage par mètre linéaire, et sx > 0 et sy > 
0 les contraintes de traction dans les lits d’armatures. 

Ainsi, l’effort résistant total des armatures, pour la direction Q, est égal à : 
R(Q) = AXsxcos?Q + AYsysin?Q (5) 


L’équation (5) est valide pour les lits inférieurs ou supérieurs, avec les efforts 
résistants Rb(Q) ou Rt(Q). La condition d’équilibre pour les lits inférieurs et 
supérieurs s’écrit donc : 
AXsxcos?2Q + AYsysin?2Q > F(Q) (6) 
Chaque inéquation pour un angle Q définit une droite dans le plan (AX, AY). 
L’ensemble des lignes définit le domaine de validité dans lequel AX et AY 
peuvent être choisis. Les contraintes sx et sy peuvent varier avec langle Q . 
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Domaine de validité prenant en compte tous les angles Q 
La solution optimale est obtenue en minimisant la section totale AX + AY en 
tenant compte des contraintes imposées par l’équation (6) et AX > 0 et AY > 0. 
Cet optimum, situé sur une intersection entre deux lignes, est calculé par le 
logiciel STFER [5]. 
Le post-traitement effectué couvre l’ensemble des pas de temps du calcul. Les 
résultats enveloppes sont présentés sous forme de cartographies donnant les 
sections d’armatures requises. La Figure 5 illustre le formalisme des résultats du 
dimensionnement du des armatures. Elle présente l’exemple de la cartographie 
des sections optimales des lits d’armatures du côté extérieur du voile Est, dans 
les directions horizontale (AX) et verticale (AY). Les sections d’armatures sont 
indiquées en cm?/m. 
La cartographie des lits horizontaux met bien en évidence la nécessité de 
renforcer particulièrement la base des voiles entre les poteaux. Seuls les pieds de 
poteaux sont solidaires des pieux. La flèche de la partie inférieure des voiles doit 
donc être reprise par un ferraillage plus important que sur les autres zones, voire 
par des longrines entre les pieds de poteaux. 
Les voiles sont encastrés en tête à la dalle de toiture. La reprise du moment de 
flexion à la jonction entre dalle et voiles nécessite un renforcement plus 
important par les lits d’aciers verticaux au voisinage de la tête de voile. Il en est 
de même à la jonction entre le voile et le toit du petit bâtiment chicane, sur 
lequel le bâtiment principal prend localement appui sous l’effet de la 
surpression. 
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Cartographie des sections d’armatures extérieures du voile Est 
Pour le calcul des armatures transversales, les deux composantes d’effort 
tranchant permettent de définir un effort tranchant pour chaque facette dont la 
normale fait un angle 0 par rapport à l’axe local x : 

QQ = QxcosQ + QysinQ (7) 

L’effort tranchant maximal ainsi obtenu est utilisé pour la justification avec les 
formules de l’Eurocode 2 

e Optimisation par une méthode non linéaire : 
Le ferraillage est dimensionné jusqu'ici par une approche linéaire. Or la réalité 
du comportement du béton armé fait intervenir un certain nombre de non- 
linéarités. Ceci concerne en premier lieu les lois de comportement des 
matériaux. L’acier a un comportement élastoplastique et une faible déformation 
plastique ne remet pas nécessairement en cause la stabilité de la structure. La 
résistance du bétondépend du tenseur de pression hydrostatique et du tenseur des 






































contraintes déviatoriques. Le comportement du béton est complexe, variant d’un 
caractère fragile en traction peu confinée à ductile en compression très confinée 
[7]. Enfin, l’interaction entre le béton et les armatures d’acier constitue une autre 
non linéarité dans le comportement du béton armé. La prise en compte de ces 
non-linéarités dans le dimensionnement permet d’adopter des critères de 
dimensionnement moins conservateurs, et donc de diminuer les sections 
d’armatures. 
Pour traiter ces non-linéarités avec LS-DYNA, nous avions imaginé initialement 
de mettre en oeuvre un modèle plaques multi-feuillets et poutres multi-fibres [3]. 
Les coques multi-feuillets permettent d’attribuer une loi de comportement sur 
chaque « rangée » de points d’intégration dans l’épaisseur de la coque. Le 
ferraillage est alors modélisé comme une plaque de section équivalente à celle 
d’un lit d’armatures. Pour les poutres multi-fibres, la discrétisation de la section 
en fibres qui permettent ensuite de modéliser les armatures nécessite un 
développement externe. Bien que cette méthode soit avantageuse en termes de 
temps de calcul et de quantité de mémoire requise, elle présente des limites du 
point de vue de l’objectif d’optimisation. Dans la pratique « d’Ingénieur », 1l est 
relativement complexe de modéliser des poutres béton armé en retombée de 
dalle. La réglementation préconise la mise en oeuvre de poutres dont la section 
est un « T » avec des longueurs d’ailes (parties de la dalle) forfaitaires. Le 
questionnement « d’Ingénieur » porte sur les parties d’efforts reprises d’une part 
par la poutre proprement dite et d’autre part par la dalle. Dans le cas d’un calcul 
hautement non linéaire, cette répartition des efforts est très complexe et ne peut 
être représentée pour un modèle coque avec des parties en retombée. 

C’est pourquoi l’optimisation par calcul non linéaire est traitée en modélisant 

explicitement le béton par des éléments volumiques (hexaédriques) et les 
armatures par des éléments poutres (Figure 
7). Un couplage cinématique est défini entre 
les noeuds des éléments de béton et les 
noeuds des éléments d’acier via la fonction 
*CONSTRAINED_LAGRANGE_IN_SOLI 
D [8]. La loi *MAT_CSCM_CONCRETE 
est utilisée pour modéliser le comportement 
élasto-plastique endommageable du béton, 
car elle permet de tenir compte d’un grand 
nombre de possibilités de sollicitations du 
béton. [7]. Le comportement de l’acier est 
modélisé par la loi 

*MAT_PLASTIC_KINEMATIC, car elle 
permet de calculer les déformations plastiques des armatures. 
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Comme nous l’avons expliqué auparavant, ce type de modélisation exclut les 
phénomènes de torsion. L’optimisation des sections de ferraillage porte donc 
uniquement sur les armatures de flexion. Le critère de vérification est extrait du 
règlement UFC 3-340-2 qui préconise une rotation maximale aux appuis de 2 
degrés [9]. 

Cette démarche d’optimisation permet d’obtenir selon les zones un abattement 
des quantités de ferraillage flexion de 10% à 70% par rapport à la méthode 
linéaire. Le ferraillage à mettre en oeuvre doit cependant respecter le minimum 
forfaitaire requis par l'UFC 3-340-2. La Figure 10 montre ainsi les ferraillages 
des voiles, dalles, poteaux et poutres obtenus avec la seule optimisation (à 
gauche) et les ferraillages enveloppes après la prise en compte des minima 
UFC (à droite). Les rotations aux appuis calculées avec le ferraillage optimisé 
n’excèdent pas 0,3 degrés, mais des sections inférieures induisent une rupture 


des armatures. 
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V. Des perspectives : 


Pour dimensionner les armatures d’un bâtiment en béton armé à une explosion 
interne, nous avons opté pour une méthode de calcul linéaire couplant un calcul 
aux éléments finis en dynamique et les règles de l’Eurocode 2. Les calculs aux 
éléments finis ont été réalisés à l’aide du logiciel LS-DYNA, en modélisant les 
voiles et les dalles par des éléments coques et les renforts par des éléments 
poutres. Les ferraillages des dalles et des voiles ont ensuite été calculés à l’aide 
du logiciel STFERO. Les ferraillages des poteaux et des poutres ont été établis 
par la méthode du diagramme d’interaction. 
Les avantages de ce type d’approche par rapport à l’utilisation d’une approche 
statique équivalente sont multiples. La variation temporelle et spatiale du 
chargement a pu être prise en compte de manière réaliste. Les ondes de 
pression caractérisant les effets de souffle se recombinent entre elles au fur et à 
mesure de leurs réflexions sur les parois. Elles présentent donc une multitude 
de phases positives et négatives, d’amplitudes et de durées variées. La 
modélisation de la totalité du bâtiment a permis de déceler les phénomènes 
physiques résultant des singularités des géométries et des profils de pression. 
Ainsi, une analyse de la dépendance des phénomènes de torsion aux conditions 
de liaison des pieds de poteaux aux têtes de pieux a permis de définir les 
dispositions constructives à mettre en oeuvre. Sur l’ensemble du bâtiment, cette 
méthode a permis de limiter les zones de ferraillages denses. 
Cependant, l’un des avantages majeurs des méthodes de calculs avancées est la 
prise en compte des non-linéarités du comportement des matériaux, qui permet 
de retenir des critères de vérification des structures moins conservateurs. La 
seule approche linéaire ne permet donc pas d’obtenir un ferraillage optimisé. 
Une approche non-linéaire a été proposée, avec la modélisation des non- 
linéarités des comportements du béton armé. Cette approche nécessite toutefois 
des ressources bien plus importantes car le béton est modélisé en éléments 
volumiques et les armatures sont explicitement modélisées en éléments poutres, 
avec un couplage de leurs comportements. L’optimisation du ferraillage a donc 
été restreinte à une section courante du bâtiment et n’a concerné que les 
armatures de flexion. La modélisation partielle du bâtiment prenant en compte 
des hypothèses de symétrie ne permet pas de restituer tous les modes de réponse 
du bâtiment, en particulier les effets de torsion. Le critère de vérification 
préconisé par UFC 3-340-2, qui porte sur les rotations aux appuis, est un 
critère global. Le passage à un dimensionnement de la structure sur un critère 
local, c’est-à-dire un critère qui soit considéré à l’échelle du matériau et non plus 
à celle de la structure, apporterait une plus-value supplémentaire à l’utilisation 
d’approches non linéaires. L’optimisation serait ainsi réalisée en intégrant des 
critères de rupture du béton et des armatures. Malheureusement, 1 UFC 3-340-2 
ne fournit pas ce type de critères locaux mais seulement des critères globaux de 
rotation de section. 






































Sur ce type de dimensionnement par des méthodes de calculs avancées, pour 
une approche linéaire ou non-linéaire, la modélisation de la réponse du 
bâtiment peut être améliorée par la prise en compte des raideurs du sol au 
niveau des têtes de pieux, à conditions que celles-ci soient caractérisées au 
préalable. L’assouplissement des conditions aux limites aux pieds de poteaux 
peut en effet influer localement sur les dispositions constructives à mettre en 
œuvre. 


VI. Application au Maroc : | 
Dans son article intitulé CONSTRUCTION DES BATIMENTS 


OPTIMSATION DES COUTS ET DES DELAIS Abdelhak EL MOULI 
ingénieur d’état en génie civil, lauréat de l’Ecole Mohammadia des Ingénieurs, a 


abordé le problème de ferraillage :  L’une des principales raisons pour 
lesquelles le façonnage du ferraillage génère des quantités importantes de chutes 
est le fait que les ferrailleurs s’efforcent systématiquement de couper les barres 
longitudinales des élévations afin d’avoir des attentes au niveau de chaque 
plancher, alors qu’il est tout à fait possible de les conserver sur toutes leurs 
longueurs ou du moins réduire le nombre de coupes. Quoique cette méthode 
nécessite un système de stabilisation pour équilibrer le ferraillage, elle permet de 
faire de grandes économies. D'ailleurs, il est recommandé de toujours privilégier 
le recours aux chutes plutôt que l’utilisation de barres saines. Quant aux 
élévations, le respect strict des plans de ferraillage contribue à la diminution du 
gaspillage du matériau. En revanche, la variation des espacements des armatures 
transversales est de nature à créer la confusion chez les ferrailleurs et donc des 
erreurs, raison pour laquelle il est souvent préférable de les harmoniser en 
optant, évidemment, pour les plus petits. Cette petite adaptation, aussi simple 
soit-elle, est en mesure de soulager grandement les équipes et les aider à 
accélérer le rythme de l’assemblage. 


VII. Réalisation du projet d’optimisation : 


A mon avis , la réalisation de ce genre de projet d’optimisation lié au 
ferraillage sera plus efficace si il est effectué par un bureau d’études puisque 
cela sera moins cher en terme de budget une grande société va demander un 
cout élevé , d’autre part c’est le bureau d’étude qui s’occupe de tous les 
étapes de réalisation d’un ouvrage y compris les plans béton armé et le 
ferraillage. 
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